Fluorescence Intensity Multiple Distributions Analysis: 
Concurrent Determination of Diffusion Times and Molecular Brightness 



The present invention relates to a method for characterizing samples having fluores- 
cent particles, said method being based on the concurrent determination of diffusion limes 
and specific brightnesses, or measures thereof. 

Fluorescence correlation spectroscopy (FCS) has proven to be a technique with single 
molecule sensitivity. Magde et al, (Phys. Rev. Lett. 29: 704 - 708, 1972) demonstrated the 
feasibility of detecting molecular number fluctuations by fluorescence correlation spectros- 
copy- Since then an increasing number of publications has appeared, aimed at improving 
the performance and accuracy of this technique. A progress was the implementation of con- 
focal detection optics (Koppel et aL, Biophys. J. 16: 1315 - 1329, 1976; Rigler and Widen- 
gren, Bioscience 40: 180 - 183, 1990) and the use of silicon photon detectors (Rigler et aL, 
Eur- Biophys. I 22: 169 - 175, 1993). This development pushed the detection limit below 
the single-molecule level (Rigler et aL, In: Fluorescence Spectroscopy: New Methods and 
Applications. O, S. Wolfbeis, editor. Springer.. Berlin. 13 - 24, 1993; Eigen and Rigler, 
Proc. Natl. Acad. Sci. USA h : 5740 - 5747, 1 994; Brand et aL, J. Phys. Chera. 101 : 43 13 
-4321, 1997; Eggeling eta]., Proc. Natl. Acad. Sci. USA 95: 1556- 1561, 1998). 

In recent publications comprising fluorescence fluctuation spectroscopy the attention 
has been drawn towards analyzing the histogram of the number of photon counts rather 
than the autocorrelation function (Qian and Elson, Biophys. I 57: 375 - 380, 1990; Fries et 
aL, J. Phys. Chem. 102: 6601 - 6613, 1998; Chen et aL, Biophys. J. 77: 553 - 567, 1999; 
Kask et al., Proc, Natl. Acad. Sci. USA 96: 13756 - 13761, 1999).Lately f FCS has found a 
complement in the form of fluorescence intensity distribution analysis (FIDA). FIDA has 
been introduced as a method for analyzing mixtures of fluorescent particles. It is based on 
the detection of instantaneous photon emission rates from an open confocal volume. Cen- 
tral part of the method is the collection of photon count numbers, recorded in time intervals 
of fixed duration (time windows) and using this information to build up a count number 
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' histogram, A theoretical probability distribution of photon count numbers is fitted against 
the obtained histogram yielding specific brightness volucs, q s as well as concentrations, c, 
Tor all different species in ihe sample. The term "specific brightness" generally denotes ihc 
mean count rate of the detector from light emitted by a panicle of given species situated in 
a certain point in the sample, conventionally in the point where ihc value of the brightness 
profile function is unity. 

The historic predecessor of FIDA is FCS, which distinguishes different species on 
die basis of their characteristic diffusion times t, by analyzing the second order autocorre- 
lation function of light intensity, G(t) « (/(0)/(0)-{/) 2 . Parameters which can be deter- 
mined by FCS (in addition to diffusion times x) are not, however, concentrations and spe- 
cific brightness values of different species separately, but products of the form cq 1 . While 
fluorescence intensity distribution analysis (FIDA) relies on a collection of instantaneous 
values of the fluctuating intensity, FCS analyzes the temporal characteristics of the fluc- 
tuations. Hence, the two methods represent complementary tools: FCS resolves compo- 
nents with different diffusion coefficients, while FIDA distinguishes the species according 
to their different values of specific molecular brightness. 

The present invention is aimed at developing a method that not only unifies the ad- 
vantages of both techniques, but rather yields a synergistic effect. 

According to the present invention, a new fluorescence fluctuation method is intro- 
duced that is based on the global analysis of a set of photon count number histograms, re- 
corded with preferably multiple widths of counting time intervals simultaneously. This 
Fluorescence Intensity Multiple Distributions Analysis (FIMDA) distinguishes fluorescent 
species on the basis of both, the specific brightness and the translational diffusion time. 
The combined information, extracted from a single measurement, increases the read-out 
effectively by one dimension and thus breaks the individual limits of prior art techniques, 
i.e. FCS and FIDA. FIMDA can be widely applied for monitoring molecular interactions 
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including receptors and tigands or antibodies and antigens, which arc both of great rele- 
vance in the life sciences. 

The present invention describes the dependence of photon count number distributions 
on diffusion coefficients and applies it to a series of photon count number histograms cor- 
responding preferably to different widths of counting time intervals. In the following, 
sometimes the wording "experimentally determined distribution function'* or "experimen- 
tally determined distribution" will be used instead of "histogram". As explained, the key is 
to analyze a set of distributions which is sensitive to the translalional diffusion of particles. 
FCS detects the dynamics of particles, because it compares the instantaneous intensities at 
time intervals separated by a certain delay. In order to make the distribution of photon 
count numbers sensitive to the temporal evolution of intensity one may alternatively - ac- 
cording to the present invention - choose to build a set of photon count number distribu- 
tions corresponding to different time windows. The choice of the time windows should 
preferably span a range comparable to the delay values used in FCS. 

In contrast to other two-dimensional FIDA techniques (Kask et al., Biophys. J. 78, 
2000), which use two detectors, in FIMDA only a single detector is needed, though multi- 
ple detectors might be used. 

While the ability of the method according to the present invention to determine spe- 
cific brightness values 3 diffusion times ) and concentrations from mixtures is demonstrated 
on simulated data, its experimental utilization is shown by the determination of the binding 
constant of a protein-ligand interaction exemplifying its broad applicability in the life sci- 
ences. Since the values for all three parameters, i.e. brightness, diffusion and concentration, 
can be attributed to each species, the method according to the present invention is a more 
powerful tool of analysis than applying FCS and FIDA separately. The statistical accuracy 
of each extracted parameter determined according to the present invention is comparable to 
that in FCS or FIDA, often even better. 



In the following, a modification of the theory of FIDA is presented, which is a suit- 
able approximation for experimental purposes. In FIDA, a convenient representalion of a 
photon count number distribution P(n) is its yencratinjj function, defined as 



(1) 

The simple theory of I : 1DA assumes (i) thai: molecules are immobile during the 
counting time interval, and (ii) that the light flux from a molecule can be expressed as a 
product of a spatial brightness function B(r) (this is a function of spatial coordinates of the 
molecule characterizing the equipment) and a specific brightness q (characterizing a certain 
molecule species). Under these two assumptions, the distribution of the number of photon 
counts, emitted by molecules from a volume element dV, is double Poissonian and the cor- 
responding generating function reads 



where £ is the complex argument of the generating function, c is the concentration of mole- 
cules, and T is the width of the counting time interval. The representation is particularly 
convenient, because contributions from independent sources, like different volume ele- 
ments or species, are combined by simple multiplication of the contributing generating 
functions. The generating function of P (n) for a single species is 

^,(0 = exp[cJ(e«-^(r)r M )^ ) ( 
while accounting for multiple species simply yields 



(4) 
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Ihc integral on the righl hand side of Eq. 4 is calculated numerically, bui instead of 
the three-dimensional integration over spatial coordinates, a one-dimensional integration 
coordinate x = ln|* fl /JS(r)] is introduced. The relationship between the brightness B and 
Ihc coordinate x is therefore B(x) = V In FIDA il is- suitable to express the function 

(Wldx, which describes the brightness profile in one-dimensional representation, by (he 
formula: 



dV 

— = A 0 (l + a { x + a 1 x 1 )x°> (5) 



Here a,, a 2 and a 3 arc empirical adjustment parameters granting for a sufficient flexi- 
bility lo fit the measured histograms with high precision. However, a representation in the 

f dV _ A / 2 3 

X 0 ( * + a,X + " lX > mi « ht ^ bc us6d - ^ se > e clion of coefficient 4, and B a 

is nothing but the selection of the units of P and B. Usually, they are determined from the 
conditions 



\BdV = \, 

(6) 
(7) 

So far, a simple version of the theory of FIDA has been described. For the purposes 
of FIMDA according to the present invention, one has to abandon the assumption that 
molecules are immobile during die counting interval. Surprisingly, according to the present 
invention one does nor abandon Eq. 2, as well as the following equations, but o.ne will re- 
define the meaning of some variables instead, x is still a variable related to the spatial 
brightness profile, but now it characterizes the path of the molecule rather than its position. 
B is the spatial brightness averaged over the path rather than determined at a fixed position 
of the molecule. V is not the volume in space but dV/dx still expresses the probability that a 
molecule has a given value of x. If one would keep the original meaning of c and q. one 
would have to d velop a theory predicting how^ 0 , fl[ , a , and a 3 depend on the counting 



'lime iiucrval 7". However, according lo the present invention another approach has been 
chosen. The normalization conditions have been kept (Eqs. 6 and 7) and it was even found 
possible to apply a single selection of ihe values /f 0 , a u <v 2 and n 2 for n set at different time 
windows. The consequence of this selection is that in liqs. 2-4 c is an apparent concentra- 
tion (c app ) and q is an apparent brightness (q m ) which both depend on the width of the 
counting time interval T. 

In ihe following, a theory is presented predicting how and q m depend on T. The 
case of single species is studied and the first and the second factorial cumulants of the dis- 
tribution corresponding to Eq. 3 are calculated. The factorial cumulanls are defined as 



K *={i]w®i<-> (8) 



yielding: 



*.=(»)=.^^,7 (9) 

where normalization conditions given by Eqs. 6 and 7 have been used. (Note that Eqs. 9 
and 10 are in total agreement with Qian and Elson's formulae (Biophys. J. 57: 375 - 380, 
1 990) derived under the assumptions i and ii.) From Eq. 9 one can conclude that 

c «*V)<la PP {T) = (l), (11) 

where (/) ■ (n) T IT is the mean count rate, which does not depend on the choice of T. One 
shall proceed by employing the following relationship between the second cumulant of the 
count number distribution P(n;T) and the autocorrelation function of fluorescence intensity 
G(0 = (/(0)/(0)-(/)\ 
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Introducing (he nolaiion 



r(r,= ^^F/ rf ''. HG( ' , "'' ) - (,3) 

one gets from Eqs. 12 and 10 

<w 0X, 2 (7) - e m (0)^ 2 (0)rCH. (14) 
From Eqs. ] I and 14 one gets 

**.(7)-«,,<o)r<r), (16) 

As the concluding step in the theory underlying the present invention, one shall sub- 
stitute the expressions of G(t) from FCS into Eq. 13. If one ignores triplet trapping and 
studies pure diffusion, then e m (0) is the true concentration c , and q m (0) is the true spe- 
cific brightness q . Applying a Gaussian brightness function (Arag6n and Pecora, J. Chem. 
Phys. 64: 1791 - 1803, 1976), the autocorrelation function is 



(17) 



denoting /; as the diffusion coefficient and a,- as ihc radial and a z as the longitudinal dis- 
tance, where the Gaussian profile has dropped e m times. The integrals in Uq. 13 yield the 
correction factor for tran.$lational diffusion 



/J (J +r)ananh 



(18) 



where r = cr r 2 and /3 = a/ /a / . For reasons explained below it is useful to calculate 
the first order terms in Eq. I &: 



,dt(i i 

l+ -r\~T + — 
6 \ a > 9 * Jj 



-i-i 



+ 0(D 1 ). 



(19) 



However, from theoretical considerations as. well as from measurements it is known 
that simple physical models like Gaussian or else Gaussian-Lorentzian do not exactly rep- 
resent the actual brightness profile (Kask et al., Proc. Natl. Acad. Sci. USA 96: 13756 - 
13761, 1999). Therefore, according to the present invention a modified Eq. 19 has been 
used and a fitting parameter a has been introduced, that preserves the first order terms in 
Eq. 19: 



1 + 



DT 
6a 



(20) 



By matching the second order terms the Gaussian brightness profile would corre- 
spond to a - lll t but we rather choose a to be an empirical parameter, which can be de- 
termined by the fitting procedure. From Eqs. .15 and 16 one can express the apparent pa- 
rameters of a pure diffusion process: 
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(21) 
(22) 



Another phenomenon involved is thai of inicnsily fluciuaiions due lo trapping of 
molecules into a triplet excited state (Widengreij et al., J. Phys. Chem. 99, 1336b' - 13379, 
1995). To obtain a good fit, particularly at values of T comparable to the triplet lifetime 
(which is typically 2 us), an additional factor has been introduced into G(t): 



1 + Ktexp 



X 



(1 + KT) 2 



(23) 



where k is the singlet to triplet transition rate and x is the triplet lifetime. As the following 
step, one may consider Eq. 23 with an additional factor of cq 2 as a correlation function of 
an ensemble of immobile particles undergoing triplet transitions: 



(24) 



From Eq. 13 and 24 one can compute: 



\2-KZ 



l + jcr — 
T 



1-e 
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+ (1+KT) 1 



(1 + /CT) 3 



(25) 



Unlike diffusion, that influences only higher cumulants of photon count numbers, 
triplet corrections also shift the first cumulant by the factor 1 /(l + kt). 
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1 + fCT 



(26) 



Solving those equation's with i aspect to (j^md c ( t ^ ] yields 



= ' . . (27) 



C°CO-*r^or). (28) 

Now, having solved the problems with diffusion and triplet transitions separately, 
one shall study the joint problem. Usually, the time scale of triplet transitions is much 
shorter than that of diffusion. Therefore, one is justified to replace c and q in Eqs. 21 and 
22 by and q<fp . This lets one combine Eqs. 21, 22, 27, and 28 to express c m and 



c (T) = 

U " A r^cor^oxi+icr) 1 ^ 
9*(T)=«r l ^(7 , )r (V (r). 



After having derived these expressions for <r Dpp and q m , the data simulations and the 
experiments should verify their validity. 

Taking into consideration the foregoing comments, the present invention relates in a 
general aspect to a method for characterizing samples having fluorescent particles, com- 
prising the steps of: 

- monitoring intensity fluctuations of fluorescence emitted by the particles in ac least one 
measurem nt volume by detecting sequences of photon counts by at least one photon 
detector, 



U. I'M, 1 I • U J 



yuan n 11:11 1 vuh a&c 1 jbii&t iwumi 



11K UO I I 0. 



com- 



11 



- dctcrjuijiing Iron; (he sequences of phoion counts intermediate statistical daia 
prising least iwo probability functions, ^(«,). A (»,),- , oC ±tt number of photon 
counts. n,ji lt ..., detected in different sels of counting lime intervals, 

- determining from said intermediate statistical data a distribution of panicles as a func- 
tion of at least two arguments, wherein one argument is a (measure of) specific bright- 
ness of the particles and another argument is a (measure of) diffusion coefficient of the 
pailicles. 

In a preferred embodiment, the diffusion coefficient of a particle can be determined. 
It might however be also advantageous 10 use a measure which describes the diffusion of 
pailicles, e.g. diffusion time. The equivalent considerations apply to specific brightness, or 
a measure thereof. 

Generally speaking, the present invention therefore comprises the steps of 

(i) Optically sensing signals resulting from fluorescence emitted from fluorescent par- 
ticles like c.g. molecules present in the sample, and 

(ii) Electronically processing the optically sensed signals, wherein said processing 
comprises the step of concurrent determination of specific brightness and diffusion 
times of said particles based on probability functions of photon counts detected in 
different sets of counting time intervals. 

In a preferred embodiment, each set of counting time intervals consists of intervals of 
equal width while different probability functions ^fc).^ (»,).... correspond-to counting 
time intervals of different widths T {1 T 2 ,.... In principle, at least two distributions with dif- 
ferent widths of the counting time interval are needed in order to obtain information about 
the diffusion, but it might be preferred to use more than two distributions for analysis. 



In a further preferred embodiment, in each set of counting time intervals these 



inter- 
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vals arc consecutive in time. However, in each sel of counting time intervals these might 
also overlap. 

Jl is particularly preferred that said distribution function ofparticlcs is determined by 
fitting the experimentally determined probability functions ^(n,),^(n 2 ).... by 
corresponding theoretical probability functions 1\ (a, ),/> (iij),... 

In a preferred embodiment, said intermediate statistical data are processed applying 
inverse transformation with regularization (ITR) or inverse transformation with constraints 
(ITC) or inverse transformation with both regularization and constraints (ITRC). Because 
of statistical errors and limited size of measured data, inverse transformation as such is 
olien an ill-posed mathematical problem, characterized by wild oscillations in its outcome. 
ITR, ITC and ITRC stabilize the mathematical problem by looking for a "regular" (e.g. a 
smooth) or constrained solution, for example by minimizing the sum of squared deviations 
of statistical data and a function of the solution itself, penalizing "irregular", usually irre- 
producible structures in the outcome, or values having no physical meaning. An example of 
constraining is disallowing negative values for concentrations (For the method of ITR, see 
e.g. W.H. Press et al. Numerical recipes in C: the art of scientific computing, 2 nd edition, 
Cambridge University Press, 1992, p, 808, which is hereby incorporated by reference). 

In a preferred embodiment, the theoretical distributions P l (o i ),J' 1 (a 2 ) ) ... are 
calculated through their generating functions G Pla) (l) = £| B .P(n). In particular, said 

rt 

distribution function of particles can be determined by fitting the experimentally 
determined probability functions P t (n l ),P 1 («,),... by corresponding theoretical probability 

functions P t («,),;>(/»,) which are preferably calculated through their generating 

functions G, (n) (a = 2f"A«) • 



In calculations of the theoretical distributions / > l (n,)./* a (n l ),... the optical spatial 
brightness function B(r) might be accounted for by a separately determined relationship 
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between brightness B and volume clcmenis dV. Depending on the position of a panicle in 
the measurement volume, the detected values of its specific brightness might vary. Espe- 
cially in cases of high polydispcrsity of the sample aiid/or low specific brighlncss differ- 
ences of the particles in (he sample, it is preferred - in respccl i 0 the accuracy of the meas- 
ured specific brightness values - lo take the spatial brightness function of the optical set-up 
into account. The relationship between the spatial brightness B and volume elements dV 
can be expressed through a variable *=ln(V*) by a relationship 
dV _ 

— - 4,(l + fl,* + fli* ) x °>, where B 0 is maximum brightness and A 0 , a,, o 2 and q 3 arc 
empirical parameters of the spatial brightness function. 

In a preferred embodiment, the generating function is calculated using the expression 
<?«) = «P[.J*efo) fd 3 r(e<W" - 1)] , where c(q) is the density of particles with spe- 
cific brightness g, T is the length of the counting time interval, and B(r) is the spatial 
brightness profile as a function of coordinates. 

In a further preferred embodiment, concentrations of particles are selected to be ap- 
proximately one or less molecules per measurement volume. Experiments if performed at 
significantly lower concentrations than one particle per measurement volume would result 
in a slow acquisition of meaningful information because most of the data collection time is 
spent on waiting, i.e. with no particles in the measurement volume. In the other limit of an 
increasing high number of particles, photon count number distributions approach a Gaus- 
sian limit; therefore it is increasingly difficult to resolve contributions from different spe- 
cies. 

Preferably, either an avalanche photodiode or a photomultiplicr is used as a photon- 
detector. It should be remarked that utilization of a single detector to collect distributions 
using a set of different widths of the counting time interval is favourite due to its fastness. 
However, in some instances it might be preferred to use at least two photon detectors, 
monitoring fluorescence of different wavelenghts or polarization. 
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The method according lo the present invention can be preferably applied in the. con- 
ductance of homogeneous assays, i.e. mix-and-measurc assays with no need of separation 
slops. 



It can be eg. applied in diagnostics, high throughput drug screening, biomcd 
optimization of properties of molecules and identification of specific cell populations. 



lClilC, 



According to the present invention, confocal techniques arc particularly suited to 
monitor fluctuating intensity of fluorescence. They may be applied - as outlined above - to 
a wide field of applications, such as biomedicine, etc. The conjugate focal (confocal) tech- 
nique is based on using a point source of light sharply focused to a diffraction-limited spot 
on the sample. The emitted light is viewed through a spatial filter (pinhole) that isolates the 
viewing area to that exactly coincident with the illuminating spot. Thus, the illumination 
and detection apertures are optically conjugated with each other. Light originating from 
focal planes other than that of the objective lens is rejected, which effectively provides a 
very small depth of field. Therefore, in a particular preferred embodiment of the present 
invention, a confocal microscope is used for monitoring the intensity of fluorescence. In 
order to achieve a high signal-to-noise-ratio, it is useful to monitor fluorescence using an 
apparatus that comprises: a radiation source (12) for providing excitation radiation (14), an 
objective (22) for focussing the excitation radiation (14) into a measurement volume (26), a 
detector (42) for delecting emission radiation (30) that stems from the measurement vol- 
ume (26), and an opaque means (44) positioned in the pathway (32) of the emission radia- 
tion (30) or excitation radiation (14) for erasing the central part of the emission radiation 
(30) or excitation radiation (14). It might be particularly preferred to use an optical set-up 
described in detail in Figure 7. 

The advantage of FIMDA (present invention) over FCS (prior art) is that the method 
yields genuine concentrations of components in the sample, instead of the products of con- 
centration and brightness squared in FCS. Only the independent determination of at least 
one of the two specific brightness values enables FCS to determine two concentrations un- 
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ambiguously. However, inexperienced users of FCS often silently assume equal specific 
brightness when resolving two components. This assumption can cause significantly biased 
results. FIMDA brings this issue to the focus of analysis. 



Another advantage of the presented method is its versatility. If FCS (prior art) or 
FIDA (prior art) fail to detect a particular readout upon a biochemical reaction, FIMDA 
Cpresent invention) might be able to succeed. The biochemical reaction is not necessarily 
limited to the binding of two components, but can be any chemical reaction of interest. The 
possibility of using only one detector for recording two physical characteristics in a single 
measurement makes FIMDA a very efficient method of analysis which saves precious as- 
say development time. 



Because two read-out parameters are recorded in one measurement preferably with 
one detector, the present invention is an efficient and reliable method of analysis. The in- 
troduction of the generation function facilitates data evaluation and makes the method a 
fast analytical tool. FIMDA is more universal applicable to biochemical assays than FCS 
and FIDA applied separately, because each assay component can be distinguished by both 
diffusion time and brightness. The statistical accuracy of determinmg the diffusion time in 
FIMDA is better or at least as good as in FCS. 

The invention is illustrated by the following figures and examples, which are not in- 
tended to limit the scope of the invention. 

Figure I illustrates the count number distributions and fits of a 3.8 nM Cy5 solution 
recorded simultaneously at different time windows T. The weighted residuals for the dif- 
ferent time windows are shown in the lower part of the figure. 

Figure 2 illustrates the fitting results of simulated data for a mixture of 3 compo- 
nents. The simulated brightness (in kHz) and diffusion rime (in us) values for the compo- 
nents are: (30 kHz, 192 us); (120 kHz, 192 us); (120 kHz, 64 us). The contributions to the 
total intensity are 10.8 kHz, 20.4 kHz, and 14.4 kHz, respectively. The graph presents the 
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results of FIMOA from 20 independent realizations of simulations, each eoaospondiug to 
an experiment of 60 s duration, 

Figure 3 illustrates the binding of pTyr-VaI-Asn-V a l-Lys(Cy5) to SH2. The solid 
curve results from a hyperbolic lit, yielding a binding eonslnnt of A'„ = 1.54 i O.J 4 U R 

Figure 4 illustrates the experimental set-up used in Experiment 2. Radiation emitted 
by a laser passes an OD filter and reaches a dichroic mirror which reflects the radiation 
towards an objective having its focus within the sample under study. Fluorescence emitted 
from the sample passes the objective and reaches the dichroic mirror which is transparent 
for the fluorescent emission. After passing a bandpass filter and a pinhole, the emiss.on 
reaches an avalanche photo-diode used as pan of the detector. By means of a photon 
counting unit and a computer, specific brightness and diffusion can be determined accord- 
ing to the present invention. 

Figure S shows 10 count number distributions with time windows 40, 60, 120, 200, 
400, 600, 800, 1200, 1600 and 2000 us from a 0.8 nM Cy 5 solution. 

Figure 6 shows the calculated apparent specific brightness of the dye as a function of 
counting time interval/ evaluated by FJDA. 

Reference is now made to Figure 7 which shows one embodiment of an apparatus 
adapted for use in performing the method according to the present invention. Apparatus 10 
comprises a laser 12 which serves as a light source for illuminating the sample by a bundle 
of coherent monochromatic excitation radiation 14, Excitation radiation 14 is paralleled by 
a lens 16 and reaches a dichroic mirror 20. Preferably, the angle between, the optical axes 
18 and the dichroic mirror 20 is 45'. The dichroic mirror 20 reflects the excitation radiation 
14 in direction of an objective lens 22 having its focus 24 within a sample volume 26. 
Sample volume 26 and objective lens 22 are preferably separated from each other by a 
transparent cover glass 28, e.g. by me bottom of a commercially available micro-titer plate 
which houses the sample. The sample preferably includes fluoreaondy labeled molecules 
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or other panicle,. Due to excitation by an appropriate excitation radiation 14, the molecules 
or other parties present i„ the sample cmii radialion 30. Emiss.on rad.ut.on 30 passes the 
objective lens 22 and reaches the dicluoic mioor 20 which is transparent for cl „i ssiori M . 
d«t.on 30. Thereafter, emission radiation pas.es a filter .14 and n collimator lens 36 on the 
opfcal axes 32. A pinhole 38 is situated m ,hc focus of collate, lens 36. Emission rndia- 
uo,i 30 passmg the pinhole 38 reaches a further lens 40 and, thereafter, lS delected by the 
photo-detector 42. Within the pathway 0 f emission radiation 30, in particular between di- 
chroic rmrror 20 and photo-detector 42, an opaque means 44 is provided through which a 
central part of the emission radtation 30 cannot p^. This central part of the emission ra- 
dtanon 30 stems from areas on die optical axes 32 in front of or behind the focus 24 of the 
exctanon nation 14. Only emission radiation 30 that stems from the focus 24 or its di- 
rect neighborhood passes the pinhole 38 and reaches photo-detector 42. Instead of placing 
an opaque means 44 within the pathway of emission radiation 30, the pathway of excitation 
radtatton 14 is also suitable for positioning an opaque means 44. In pabular, an opaque 
means 44 can be positioned between laser 12 and dichroie mirror 20. The use of an opaque 
means 44 in the method according to the present invention as described in deta.1 herein 
improves the signal-to-noise ratio. 
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EXPERIMKNT 1 

MATERIALS AJVD METHODS' 
Experimental set-up 

The central opto! part of*, present FIM1M experimen, is „ confocl micnKcopc 

« s ,< « ustJ „, n„„„,™,,a; c ,„ Idali0 „ S|)cdroscopy (ltoppcl rt ^ Bjoiihjs j Bi5 _ 

1329. 1976: JUglcr ol iL, Bur. Biopbys. J. 22; .09- ,75, ,993). Fo, the excitation of „„„. 
rcscence, a beam from , oo nli „„ ous wave h ^ w _ 800 mW by ^ 

filters, passes a beam expander, and Is directed to rj,c miaoscope objective (UApo/340 
40x, N.A. U5, Olympus Optical Co. Ltd, Tokyo, Japan) by a dicta™ mi™ p, uorK . 
cence „ collected by to Slmc ^ ^ ^ ^ _ ^ ^ 

filter, and is focused to a confocal pinhole, which serves * reject the out-nf.focus light 
The Itgk, which passes the pinhole, is detected b, a silicon photon county m odule 
(SfCM-AQ-m, EG*G Optoelectronics, Vaudreuil, Canada,. An electronic counter »„. 
structed a, EVOTEC as a computer plug-i. elects the TTL pulses from the detector 
contmuousl, and calculates the count number histograms for all preselected widths of time 
windows (40, 60, 120, 200, 400, 600, 800, 1200, 1600, 2000 us) in reaMime from the 
32 MB onboard buffer. By feeding the detector outputs to a correlator, FCS measurements 
ran be performed in parallel with FIMTM experiments. 

to order to satisfy the spectral needs of the various fmorophores used i„ this audy 
Afferent lasers and spectral band-pass filters were employed. For Cy 5 (Amcrshan, Pharma-" 
era Biotech. Bucfc, UK) conjugated an arrangement of a red laser diode 

(CrystaJ (w*. Be,,!., as ^ ^ , ^ ^ ^ ^ ^ ^ ^ 

of 670 nm (670DF40, Omega Optical, Brauleboro. V 0 was used. In case of TAMRA (5 
carboxytetrante^rhodamine) Ubeled molecule this was an arrangement of a frequency 
doubled Nd-YAG laser (.Grccn u^, „ ^ CA; }}2 ^ ■ 

filter. 

The effective dimensions of 0* illuminated volume were Abated indirectly, using 
°" S,M " " 0,eO ^ (TAMRA, Cv5) with known diffusion coefficients. The ^ 
correlation fuoerions of diffusion were fined to E,. 17, i... assuming , 3 ^ mM s i(>Ml 
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. Gaussian intensity profile. The exact determination of the dimensions and profile would be 
very complex, since they are affected by both, the size of the laser beam as well a* the size 
of (lie confocal pinhole. However, in most cases the knowledge of the exact dimensions is 
mil necessary. 

The local beam radius was adjusted to -0.75 um by selecting an appropriate expan- 
sion factor of the original laser beam, resulting in a mean translation^ diffusion time of 
360 us for the free dye Cy5. This diffusion can be dearly observed when rising the time 
windows from 40 us to 2 ms. As can be seen in Fig. 1, the selected count number distribu- 
tions of a 3.8 nM Cy5 solution differ considerably. However, the major differences be- 
tween the distributions are due to the varying mean count number in different time win- 
dows used. Diffusion of fluorescent molecules causes only small but significant modifica- 
tions to the shape of each disuibution. 



The levels of background count rate are determined in a separate experiment on 
bidistillcd water and amount usually to 0.5 kHz. The main contributor to this non- 
fluctuating background light intensity is Raman scattering from water. 

Data simulations 

Real samples, comprising a mixture of molecules, which express deliberately chosen 
parameters (brightness values and diffusion coefficients), are difficult to prepare. There- 
fore, some evaluations of the new method were performed using simulated data. A number 
of sets of histograms for HMD A, FIDA, and correlation functions for FCS have been 
simulated according to the following algorithm. In a closed rectangular reservoir, a given 
number of molecules is initially randomly distributed over a high number (typically 
360x360^720) of discrete spatial grid points. Each molecule is subject to consequent diffu- 
sion simulation and jumps randomly by one grid unit either in x-, y- t or z-direction with a 
frequency corresponding to a given diffusion coefficient. The "focus" is located in the 
center of the reservoir, and the brightness distribution is assumed to be Gaussian in all 
three dimensions. When calculating the brightness integral from a molecule over a given 
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set of time intervals, the molecule can be randomly trapped, a* well as reload from the 
triplet excited .slate (where it is dark). Now one can calculate an array of brightness inte- 
grals over ^ llino intevall of a ^ ^ (e g ? ^ ^ ^ 

mixture. The brightae* integrals arc then converted into photon count numbers by gener 
«.ng a random 1'oisson number with the corresponding average. This step also accounts lor 
the tnaodueed by the detector becau.se the random number generator is used not only 
for dnvuuj random motion of molecules but also for simulating random numbers of de- 
tected photon, al eiven Jighl intensitjes - llic nuJom count numbers obta . ncd ^ ^ 

quently used to eajculate histograms for FIMDA and FIDA as well as the collation func- 



tion for FCS. 



Due to the Unite size of the simulation reservoir, some distortions of the correlation 
funcuon (i.e., deviate from Eq. 17) can be expected. The distortion, are in f ac , below 
the statistical noise level. Therefore, one can consider the stations to be an adequate 
tool for estimating statistical errors of the extracted parameter, For this purpose, typically 
30 realizations of experiments with a given set of molecular parameters were simulated 
from which the standard deviations and the coefficients of variation (CV) as the rat, of 
standard deviation to mean value were calculated! 



Filling 



A sen.* of simulaneously measuKd „ r !ira „ latel .„ gloW , y ^ m 

mg a Ma^arf, aigonte,. Th. fiMns ^ „ . modes , rf fc 

dcagnad fa FBA OOdc u a,., Proc. Natl. Acad. Sci. USA 96: 13756 - ,376, ,999) 
Thecal dtaft*- „ calculated Min8 „, same n ^ ^ ^ 

•V tm each spaci* „ „ MvUoal ^ ^ ^ ^ 

a. each wi^c,, catolated according to Eq,. 29. All pa™*** M( to 
d. M rafter . ft. eoui Pme „, (la.. A «, „. ^ „ ftom ^ , _ ^ 

dually detained beforehand from separate adjustment ™ * . 

*^«udie aqjustmeat expenments on pure dye solu- 
tions. 
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Xtoc/i emiail System 

The Crb2 (Slty-phosphopepdclt interaction 

Reccm antitumor research has been focused on tyrosine kinase growth factor rcixp. 
tors (Levitzki, Eur. J. Biochcm. 226: I - 13, 1994; Alejandro et al„ Curr. Top. Microbiol. 
Immunol. 213: 167 ~ 188, 1996; Furct ct al, J. Med. Chem. 41: 3442 - 3449. 1998). A 
critical link in the signal transduction patliway of this receptor is the interaction of its phos- 
photyrosine residue (pTyr) with the Src-homology 2 (SH2) domain of the adapter protein 
Grb2 (growth factor receptor-binding protein 2). For the recognition, a minimal peptide 
sequence of the receptor (pTyr-Val-Asn) is sufficient (MiiUer et al., J. Biol. Chcm. 271: 
16500 - 16505, 1996; Gram ct al., Eur. J. Biochem, 246: 633 - 637, 1997; Furct et al., J. 
Med. Chem. 41: 3442 - 3449, 1998). The binding partner of this peptide motive, the SH2 
domain of Grb2, can fold into a functional protein module independent of neighboring se- 
quences (Booker et al., Nature, 358: 684 - 687, 1992; Overduin et al., Cell 70: 679 - 704, 
1992). Therefore, as a model system, the bare SH2 domain (14.3 kDa) has been chosen to 
interact with a fluorescently labeled phosphopeptide (pTyr-Val-Asn-VaM,ys(Cy5)) (1387 
Da). 

The SH2 domain of Qrb2 was prepared as described elsewhere (Lowenstein et al., 
Cell 70:431 - 442, 1992; Baumann et al., Eur. J. Immunol. 24; 1799 - 1807, 1994; Miiller 
et al., J Biol. Chem. 271: 16500 - 16505, 1996). The phosphopeptide was synthesized 
using manual Fmoc solid phase chemistry and labeled with Cy5-NHS via a Lysine residue. 
An additional Valine was introduced to minimize possible interactions of the dye with the 
main recognition motive pTyr. The final compound, pTyr-VaI-Asn-Val-Lys(Cy5) was 
characterized by mass spectrometry (LC/MS, and MALD&TOF), UV/VIS, and fluores- 
cence spectroscopy. 



RESULTS 

Data simulations and test experiments 

At first, a series of measurements on a 1 ruM TAMRA solution was performed col- 
lecting data in parallel- for F1MDA as well as for FCS. This series of experiments, with 
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duration of 2 s each, was repealed in simulations using similar molecular parameters. The 
purpose of these experiments was to verify whether simulations arc a reasonable model of 
real experiments, in particular whether data simulations arc u reasonable means of predict- 
ing statistical errors of estimated parameters. The coefficients of variation of the parame- 
ters extracted from simulated data indeed coincide with the results of the real experiment, 
as can be seen in Tab. 1. 

Another scries of test experiments was repeated in a significantly shorter time do- 
main with the goal to compare FIMDA and FCS in their ability to estimate parameters of 
the triplet component A set of counting time intervals of 2 3 4, 8, 16, 32, 64, 128, 256, 512, 
and 1024 |xs was selected for this purpose. The duration of these experiments was 16 s. The 
results, presented in Tab. 2, indicate that the values for the triplet parameters estimated by 
FIMDA have similar dependence on the excitation intensity to the FCS results. The 
FIMDA results are slightly biased and have higher CV values compared to FCS, since the 
estimation of triplet parameters in FIMDA is indirect, because the shortest time window (2 
jis) is equal to the triplet lifetime. However, the main purpose of the triplet correction in the 
model is not to determine the triplet parameters, but to improve the quality of the fit and to 
remove a source of bias in the. brightness and diffusion parameters. 

Out of curiosity,' histograms for FIMDA were simulated for three-component analy- 
sis. Two of the components had equal brightness values (120 kHz), and another pair had 
equal diffusion times (192 \\s). Due to the larger number of free parameters, the simulated 
duration of experiments was increased to 60 s, so that the variations of fitted parameters 
stayed in reasonable limits. In this test, all parameters were subject to fitting. The results 
are presented in Fig. 2 as vertical bars in a plane with brightness and diffusion time as x-y 
coordinates, and the ordinate displaying the contribution to the intensity, i.e,, the product of 
concentration and brightness. The three components are clearly resolved, since the scatter 
in the location of individual bars is much smaller than the distance between the groups, 
which correspond to different components. 



Note that with FIDA alone, the components with equal brightness cannot be resolved, 
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* res .,„,«, te m , h ^ diffusioo timc ^ unreM|vcj 

biochemical system 

T„c expeW,e„„ u[ili , ation of thc „ cw ^ ^ ^ ^ 

..a(,o« of the binding constant of the above introduced Glb , CSIP ) nhos , 7 
aciion Korihi«n • ^' W-phosphopcplide inter- 

vlwc " - « » PTy r -V.,. As „. 

0.0 , 0, 0 3, ,. 3, ,0, ,0, ,00. „ ,30 /iM). Al, cxpeH™,, WHe ^ ^ 

/ 8, 50 mM NaCl, and 0.05% Pluronic- T=20°C1 ,a ♦ - . 

' CJ ' 3X1(1 a data acquisition time of 30 s ner 

^ ^ " - W * >*>. m 400, 600, 800, I200 lm " ' 

(ft 3 1.7 ± 0J fcfe we,, d^mtad fto m , ^ component ^ . wlied to to „ 
soluMon. AddiU M of excess SH2 mn P 

bo. by a longeron tirac J^JX^ ^ CharaCteriZCd 

iu 8 ate. This mix^e was ^ ^ * ^1 t ^ t0 * C ° n " 

uien analyzed by all three methods (FIMDA, FIDA and Frs\ 
■ang a two-component fit with z, and/or „ fixed den, * \ > 
thissten.f, i • ^ed, deperdmg on the method. The results of 

mis step of analysis are presented in Tab 3 It M „ i. ■ 

val„«nf , 3,ItCanbcscejltha t ^1 methods yield simUar 

values of parameters for the comnlev n« ^ 

: ^=zzzz.~zr, 

exception of FIDA which h, m , * Smular ' *** 1,16 

P on or MUA, which has difficulties due to the small n(M\ hw 

brightness of the two components. : ° } ^ 01 Specific 

As the next step of the present studies, a sample with 3 uM Sf o 

~ concentration was chosen to achieve a LS«T T ^ 

mixture of approximat ly equal proper- 
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tions of complex and free conjugate. Since it is rather difficult to resolve components with 
only a twofold difference in diffusion coefficient and even smaller difference in specific 
bri,»|,i„ css , here al so (he diffusion, lime and brightness or the complex were fixed to the 
values of Tab. 3. With the molecular parameters fixed, the conccnlralions were reliably 
determined by all methods. The results of this step of analysis arc summarized in Tab. 4. 

In the same manner, the whole series of SH2 concentrations was fitted. Pig. 3 shows 
the calculated traction bound ( Ccoaiplex / ( Ccomplcx -i- CconjUfiMB )) for FIMDA with the solid curve 
resulting from a hyperbolic fit, that yielded a binding constant for die SH2-phosphopeplide 
interaction of K D - t .54 ± 0. 14 uM. Comparable binding curves were also obtained by FCS 
and FIDA (data not shown), with K D values of 2.16 ±0.19 uM and 1.60 ± 0.19 uM, re- 
spectively. 



The data of Fig. 3 demonstrate that FIMDA is a suitable method for monitoring the 
formation of a molecular complex. FCS and FIDA experiments yielded similar * 0 values 
for this particular SH2-phosphopeptide interaction. In the literature the affinity is reported 
to vary by several orders of magnitude, depending on the peptide sequence (Mailer et al., J. 
Biol. Chen,. 271: 16500 - 16505, 1996; Gram' et al., Eur. J. Biochem. 246: 633 - 637, 
1997, Furet et al., J. Med. Chem. 41 ; 3442 - 3449, 1 998). High affinities are in the range of 
Ko - 10 - 100 nM. However, with a lysine (and C y 5 attached to it) at the +4 position of the 
phosphopeptide (defining p-Thr as the 0 position with «+' continuing on the C and on 
the N terminus) the affinity decreases to the micromoJar range. This result agrees well with 
the importance of lipophilic groups attached to 'appropriate' positions on the C-terminus, 
increasing the binding constant to the SH2-domain (Furet et al., J. Med. Chem. 41 : 3442 - 
3449, 1998). For example Val (at position pPyr+3) is making van der Waals contact with a 
large hydrophobic area on the SH2-domain. 

One of the surprising results of this study according to the present invention is that in 
each of the experiments, the statistical accuracy of the diffusion time estimated by FIMDA 
is as good as or even better than that estimated by FCS. This is a counter-intuitive result 
because FCS is directly focused on fitting a diffusion-dependent correlation function G(t), 
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while in J-IMDA the diffusion lime is estimated only indirectly, namely tluough the de- 
pendence of (he apparent brightness on the width of (he lime window. 

A further observation in this respect is lhal the CV values of the diffusion times are 
in general higher than those for the brightness values. This also holds true for the theoreti- 
cal simulations 3 „d therefore reflects an effect rooting in the measuring principle. The phe- 
nomenon can be explained qualitatively by the different ways how these quantities are de- 
termined. For simplicity, one may imagine an observation volume with a constant bright- 
ness profile B(r) inside. In this case, one only needs to measure the average count rate of a 
molecule that enters the volume to determine its specific brightness. This requires the de- 
tection of many photons per given time interval but can in principle be achieved from a 
single passage. On the other hand, for estimating the diffusion time, one has to determine 
the mean duration of the diffusion driven passage, which inevitably requires averaging over 
many events, even though many photons may be detected each time. Therefore, in an ex- 
periment of fixed duration, the specific brightness of a molecule can in principle be deter- 
mined with a higher accuracy than its diffusion time. 



EXPERIMENT 1 

Figure 4 shows^the experimental set-up used in this study. As an example for a typi- 
cal measurement Figure 5 shows 10 count number distributions with time windows 40, 60, 
120, 200, 400, 600, 800, 1200, 1600, 2000 us from a 0.8 nM Cy 5 solution. Figure 6 
shows the calculated apparent specific brightness of the dye as a function of counting time 
interval (time window), evaluated by FID A. Since the diffusion time of the dye is long 
(compared to 40 us) it spends only a short time in the confocal volume and its apparent 
brightness decreases. From this information the diffusion time can be calculated from a 
global fit. 

The viability of the present invention was tested by measurements characterizing a li- 
gand-protein interaction. The minimal peptide sequ nee pTyr-Val-Asn (originating from 
tyrosine kinase growth factor receptors) binds to the SH2-domain of the adapter protein 
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• Grb2. The conjugaic P Tyr-Va.-Asn-Val-Ly S (Cy5) was synthesized a S binding partner for 
SH2. All measurements were performed under identical conditions (e.g. same buffer: ster- 
ile filler water, 50 mM Na-phosphatc buffer pH 7.8. 50 mM NaCI and 0.05 % Pluronic; T 
- 20 -C). To obtain the K D value for the ligand protein interaction a titration experiment 
was carried out. The pTy,-V a |.A S n-Val-Ly S (Cv5) concentration was kept constant ar 0.« 
nM, while SI12 was subject to thralion (0.01. 0.03, 0.1, 0.3, 1, 3, 10, 30, 100, and 130 u M ). 
Each sample was measured 30 times keeping the acquisition time at 30 seconds per meas- 
urement FCS (prior art), FIDA (prior art) and F1MDA (present invention) data were col- 
lected simultaneously in order to compare all three methods. Table 5 summarizes the re- 
sults. It becomes clear that FIMDA is able to resolve both, the specific brightness and the 
diffusion time of both species (ligand and complex). Therefore, FIMDA is a more versatile 
tool than PCS or FIDA and will save assay development time. Comparing the standard 
deviations of all methods in determining brightness and diffusion time FIMDA is better or 
at least as good as FCS or FIDA. 
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TABLES 



Tabic 1 Comparison of coefficients of variation of - estimated parameters from series 
of experimental and simulated histograms by FIMDA. and correlation functions by FCS. 



Parameter 



Mean value 



CV (%) 



FIMDA . 



FCS 









Experimental 


Simulated 


Experimental 


Simulated 








data 


data 


data 


data 




Brightness q (kHz) 


115 


2.5 


2.3 


2.6 


2.5 


4) 

" % l 


Concentration c 


0.73 


3.1 


3.8 


3.8 


4.0 


i 

"•i 

41 


(molecules per con- 
focal volume) 












ii 


Diffusion time r (us) 


287 


5.4 


4.3 


7.3 


5.2 



Table 2 Triplet parameters, estimated from a series of experiments on 1 nM TAMRA 
solution by FCS and FIMDA at two different excitation intensities. Excitation wavelength 
532 nm, duration 16 s, time windows 2, 4, 8, 16, 32, 64, 128, 256, 512, and 1024 us. 



Peak ex- 






FCS 






FIMDA 




citation 


Triplet 


CV(»/ 


'o) Triplet 


CV (%) 


Triplet 


CV(%) 


Triplet 


CV (%) 


intensity, 


lifetime 




popula- 




lifetime 




popula- 




(kW/cm 2 ) 


(MS) 




tion 




(MS) . 




tion 




118 


1.98 


3.9 


0.182 


2.7 


3.12 


7.6 


0.137 


3.1 


187 


1.75 


3.7 


0.235 


1.6 


2.59 


3.8 


0.183 


2.0 
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Tabic 3 Comparison of estimated parameters and their coefficients of variation at a 
high receptor concentration (130 (i M). A scries of 30 experiments of 30 s duration each 
was evaluated by FIMDA, FfDA and RS. Brightness (in FIMDA And Fl|)A) and diffusion 
lime (in FJMDA and FCS) of the free conjugate were independently determined and fixed 
to 3 1 .7 kHz and 407 us, respectively, in diis analysis. 



Parameter 


Method 


Mean value from 


CV (%) 


CV (%) 






experiment 




from simulations 


— _ 






^conjugate 


FIMDA 


0.132* 


43* 


82* 


(molecules per 


FIDA 


0.196* 


76* 


71* 














FCS 


0 f»S9* 


no* 

99 


120 


^complex 


FIMDA 


0.618 


9.6 


8.0 


(molecules per 


FJDA 


0.555 


26.7 


14.9 


confocal volume) 












FCS; 


0.710 


7.9 


12.6 




..FIMDA 


39.5 


22 


2.3 .. 


Qoompicx (kHz) 


FIDA 


38.4 


5.5 


3.9 




FCS 


36.4 


3.7 


3.4 




FIMDA 


0.913 


6.9 


4.6 


^complex (ms) 












FCS 


0.898 


5.4 


7.2 



out of 30 realizations, in 5 to 10 cases zero conjugate concentration was yielded by the 
fitting program (negative values are disallowed). This indicates that in this particular ex- 
ample the conjugate concentration could not be properly determined. However, for further 
data analysis only and ^complex are needed. 
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Tabic 4. Comparison of Die estimated concentrations at an intermediate receptor 
concentration (3 u.M). In addition 10 the brightness and the diffusion time of the free conju- 
gate, also the brightness and/or the diffusion time of the complex wcrn fixed here to values 
shown in Tab. 3 



Parameter 


Method 


Mean value from 
experiment 


CV (%) 


CV (%) 
from simulations 


Ccoiijugaic 


FIMDA 


0.328 


11.7 


4.9 


(molecules per 


FIDA 


■0.311 


16.1 


4.9 


confocal volume) 


FCS 


0.303 


8.4 


8.0 


c eomplex 


FIMDA 


0.437 


8.2 


3.4 


(molecules per 
confocal volume) 


FIDA 


0.455 


10.0 


3.7 




FCS 


0.467 


7.8 


8.8 
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Table 5: Comparison of I'TMHA (present invention) with prior art l r CS and HIM 
with respect to tlic determination of q. t sind K» 



[kHz] [kHz] 

FIMDA 31.7±0.32 39.6±0.9 
FCS 

FIDA 32.1±0.31 43.3±2.3 



[ms] [ms] [jiM] 

0.402±0.006 1.01±0.08 1.68± 0.27 
0.402±0.013 0.867±0.057 2.26±0.28 

1.70±0.29 



